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1. Introduction 

Adaptive flight control is a potentially promising technology that can improve aircraft sta- 
bility and maneuverability. In recent years, adaptive control has been receiving a significant 
amount of attention. In aerospace applications, adaptive control has been demonstrated in 
many flight vehicles. For example, NASA has conducted a flight test of a neural net in telli- 
gent flight control system on board a modified F-15 test aircraft ' Bosworth & Williams-Havesl 
120071) . The U.S. Air Force and Boei ng have developed a direct adaptive controller for the Joint 
Direct Attack Munitions (JDAM) ISharma et all 12006 1. The ability to accommodate system 
uncertainties and to improve fault tolerance of a flight control system is a major selling point 
of adaptive control since traditional gain-scheduling control methods are viewed as being less 
capable of handling off-nominal flight conditions outside a normal flight envelope. Nonethe- 
less, gain-scheduling control methods are robust to disturbances and unmodeled dynamics 
when an aircraft is operated as intended. 

In spite of recent advances in adaptive control research and the potential benefits of adaptive 
control systems for enhancing flight safety in adverse conditions, there are several challenges 
related to the implementation of adaptive control technologies in flight vehicles to accom- 
modate system uncertainties. These challenges include but are not limited to:_l) robustness 
in the presence of unmodeled dynamics and exogenous disturbances iRohrs et all Il985h ; 2) 
quantification of performance and stability metrics of adaptive control as related to adaptive 
gain and input signals; 3) adaptation in the presence of actuator rate and position limits; 4) 
cross-coupling between longitudinal and lateral-directional axes due to failures, damage, and 
different rates of adaptation in each axis; and 5) on-line reconfiguration and control realloca- 
tion using non-traditional control effectors such as engines with different rate limits. 

The lack of a formal certification process for adaptive control systems p oses a major hurdle 
to the implementation o f adaptive control in future aerospace systems ilacklin et all 120051: 
iNguven & Tacklinl. 1201(1) . This hurdle can be traced to the lack of well-defined performance 
and stability metrics for adaptive control that can be used for the verification and validation of 
adaptive control systems. Recent studies by a number of authors have attempted to address 
metric evaluation for adaptive control systems llAnnaswamv et ail. 120081: INguven et all 120071: 
IStepanvan et alJj2009l:lYang et al.l.l2009h . Thus, the development of verifiable metrics for adap- 
tive control will be important in order to mature adaptive control technologies in the future. 


Over the past seve ral years, various model-reference adaptive control (MRAC) methods have 


been investigated ( Cao & Hovakimyan 

2008}; Eberhart & Ward), 199^i Hovakimyan et all 2001 

Tohnson et al., 2000:; Kim & Calisd. 1997 

Lavretskv, 2009; Neuven et al., 2008; Rvsdvk & Calise 

1998|; Steinberg, 1999|). The majority of MRAC methods mav be classified as direct, indirect 


or a combination thereof. Indirect adaptive control methods are based on identification of 
unknown plant parameters and certainty-equivalence contr ol sche me s derived f rom the pa- 
rameter estimates which are assumed to be their true values llloannu & Suri.fl996ll . Parameter 
identification techniques such as recursive l east-s qu ares and n eural networks have been used 
in many indirect adaptive control methods teberhart & WardLll999h . In contrast, direct adap- 
tive control methods adjust control parameters to account for system uncertainties directly 
without identifying unknown plant parameters expl icitly. MRAC methods based on neural 
networks have been a topic of great research interest llohnson et alll200(l:lKim & Calisel.ll997l: 
iRvsdvk & Calisel.ll998ll . Feedforward neural networks are capable of approxi mating a generi c 
class of nonlinear functions on a compact domain within arbitrary tolerance dCvbenkol. 19891) . 
thus making them suitable for adaptive control applications. In particular, Rysdyk and Calise 
described a neural net direct adaptive control method_ for improv ing tracking performance 
based on a model inversion control architecture dRvsdvk & Calisei.ll998l) . This method is the 
basis for the intelligent flight control system that has been developed for the F-15 test aircraft 
by NASA. Johnson et al. introduced a pseudo-control hedging approach for dealing with 
con trol input char acteri stics such as actuator saturation, rate limit, and linear input dynam- 
ics hohnson et alJ.12000 ). Hovakimyan et al. developed an output feedback adaptive control 
to addressjssuesjvitiajaarametric i uncertainties -i an^^nmpdele^ i d^namicsJHovakimyan et 
al.. l200ll) . C ap an d Hovakimyan developed an L\ adaptive control method to address high- 
gain control llCao & HovakimvarJ. 12008 ). Nguyen developed an optimal control modification 
schemeforada^tive^controHo^mgrove^stabilit^robustnessunderfastadagtationjNguyen et 
al.. l2008i) . 

While adaptive control has been used with success in many applications, the possibility of 
high-gain control due to fast adaptation can be an issue. In certain applications, fast adapta- 
tion is needed in order to improve the tracking performance rapidly when a system is subject 
to large uncertainties such as structural damage to an aircraft that could cause large changes 
in aerodynamic characteristics. In these situations, large adaptive gains can be used for adap- 
tation in order to reduce the tracking error quickly. However, there typically exists a balance 
between stability and fast adaptation. It is well known that high-gain control or fast adapta- 
tion can result in high frequency oscillatio ns which can excite u nmodeled dynamics that could 
adversely affect stability of an MRAC law llloannu & Sunl.ll996 ). Recognizing this, some recent 
adaptive control methods hav e begun to address fast adaptation. One such method is the L\ 
adaptive control iCao & Hovakimyan! 12008 ) which uses a low-pass filter to effectively filter 
out any high frequency oscillation that may occur due to fast adaptation. Another approach 
is the optimal cont rol mod ification that can enable fast adaptation while maintaining stability 
robustness iNguven et alj,[200~8 ). 

This study investigates a hybrid ad aptive fl i ght control method as another possibility to re- 
duce the effect of high-gain control iNguven et alj.l2006l) . The hybrid adaptive control blends 
both direct and indirect adaptive control in a model inversion flight control architecture. The 
blending of b oth direct and indirect adaptive control is sometimes known as composite adap- 
tation jloannu & Sunl.fl9% ). The indirect adaptive control is used to update the model inver- 


sion controller by two parameter estimation techniques: 1) an indirect adaptive law based on 
the Lyapunov theory, and 2) a recursive least-squares indirect adaptive law. The model inver- 
sion controller generates a command signal using estimates of the unknown plant dynamics 
to reduce the model inversion error. This directly leads to a reduced tracking error. Any resid- 
ual tracking error can then be further reduced by a direct adaptive control which generates 
an augmented reference command signal based on the residual tracking error. Because the 
direct adaptive control only needs to adapt to a residual uncertainty, its adaptive gain can be 
reduced in order to improve stability robustness. Simulations of the hybrid adaptive control 
for a damaged generic transport aircraft and a pilot-in-the-loop flight simulator study show 
that the proposed method is quite effective in providing improved command tracking perfor- 
mance for a flight control system. 

2. Hybrid Adaptive Flight Control 

Consider a rate-command-attitude-hold (RCAH) inner loop flight control design. The objec- 
tive of the study is to design an adaptive law that allows an aircraft rate response to accurately 
follow a rate command. Assuming that the airspeed is regulated by the engine thrust, then 
the rate equation for an aircraft can be written as 

u> = oj* + A 6j (1) 

where w = [ p q r ] T is the inner loop angular rate vector, Au> is the uncertainty in the 
plant model which can include nonlinear effects, and uj* is the nominal plant model where 

do* = F^co + F^cr+G*S (2) 

with F*, F|, G* e R 3x3 as nominal state transition and control sensitivity matrices which are 
assumed to be known, a = [ A <p Aa A/3 ] T is the outer loop attitude vector which has 

slower dynamics than the inner loop rate dynamics, and 5 = [ A S a AS e A S, ] T is the 
actuator command vector to flight control surfaces. 



Fig. 1. Hybrid Adaptive Flight Control Architecture 


Figure 1 illustrates the proposed hybrid adaptive flight control. The control architecture com- 
prises: 1) a reference model that translates a rate command into a desired acceleration com- 
mand, 2) a proportional-integral (PI) feedback control for rate stabilization and tracking, 3) 
a model inversion controller that computes the actuator command using the desired accel- 
eration command, 4) a neural net direct adaptive control augmentation, and 5) an indirect 
adaptive control that adjusts the model inversion controller to match the actual plant dynam- 
ics. The tracking error between the reference trajectory and the aircraft state is first reduced 


by the model inversion indirect adaptation. The neural net direct adaptation then further re- 
duces the tracking error by estimating an augmented acceleration command to compensate 
for the residual tracking error. Without the model inversion indirect adaptation, the possibil- 
ity of a high-gain control can exist with only the direct adaptation in use since a large adaptive 
gain needs to be used in order to reduce the tracking error rapidly. A high-gain control may 
be undesirable since it can lead to high frequency oscillations in the adaptive signal that can 
potentially excite unmodeled dynamics such as structural modes. The proposed hybrid adap- 
tive control can improve the performance of a flight control system by incorporating a model 
inversion indirect adaptation in conjunction with a direct adaptation. 

The inner loop rate feedback control is designed to improve aircraft rate response characteris- 
tics such as the short period mode and the dutch roll mode. A second-order reference model 
is specified to provide desired handling qualities with good damping and natural frequency 
characteristics as follows: 

(s 2 + 2^,pU}pS + w 2 J <p m = CpSiat (3) 

(s 2 + 2^qCOqS + W 2 J 6,n = CqS hn (4) 

(s 2 + 2£ r a> r s + « 2 J r m = c r S rud (5) 

where cp m , 9 m , and ip m are reference bank, pitch, and heading angles; <5/ at , Sj on , and S rud are the 
lateral stick input, longitudinal stick input, and rudder pedal input; Wp, lOq, and to r are the 
natural frequencies for desired handling qualities in the roll, pitch, and yaw axes; fg, and 
are the desired damping ratios; and Cp, Cq, and c r are stick gains. 

Let p m = (p m , q m = 6 m , and r m = i p m be the reference roll, pitch, and yaw rates. Then the 
reference model can be represented as 

(Vm — KpCVm fQ J (Vm^T eSc (6) 

where u> m = [ p m q m r m ] T , K p = diag (2£ p w p ,2£ 9 w ? ,2£,.(xv), K t = diag (w*, , 

C = diag (cp,Cq,C r ) , and 5c = [ t $ Ion ^ rud ] 

A model inversion controller is computed to obtain an estimated control surface deflection 
command S to achieve a desired acceleration w d as 

6 = G- 1 (u> d - Aw - ho) (7) 

where F\, h> ar *d G are the unknown plant matrices to be estimated by an indirect adaptive 
law which updates the model inversion controller; and moreover G is ensured to be invertible 
by verifying its matrix conditioning number. 

In order for the controller to track the reference acceleration w m , the desired acceleration d> d is 
computed as 

w d = u> m + KpCVe + Ki 0) e dT - u ad 


( 8 ) 


where u> e = co m — co is defined as a rate tracking error, and u ad is a direct adaptive signal 
designed to reduce the tracking error to small bound away from zero in order to provide sta- 
bility robustness. 

Because the true plant dynamics are unknown, the model inversion controller incurs a mod- 
eling error equal to 

do - do d = do - do m - KpCVe - Kj (L> e dT + u ad (9) 

but from Eq. {7j the model inversion controller is also equal to 

cd-td d =e-(F 1 - Ff) co - (f 2 - F 2 *) a - (G - G*) S (10) 

where e = Ad; is the unknown plant model error. 


Comparing these two equations, the tracking error equation is formed as 

e = Ae + Bu ad + BFico + BF 2 cr + BG5 — Be (11) 


where e = [ fg co e dT co e j is the tracking error, F) = F) — Fj, F 2 
and 

A= [ 0 1 
[ -Ki -K v 



P 2 -F*,G = G- G*, 
( 12 ) 


The direct adaptive signal u ad is computed from a single-layer sigma-pi neural network 


Uad = W t Y (13) 

where W 6 R mx 3 is a neural network weight matrix, and Y = [ Q C 2 C 3 C 4 C 5 ] T € 
R mxl is a basis function with Q, i = 1, ... ,5, as inputs into the neural network consisting of 
control commands, sensor feedback, and bias terms; defined as follows 

C l = V 2 [ u ; T uco T pco T ] (14) 

C 2 = V 2 [ 1 a p a 2 p 2 up ] (15) 

C 3 = V 2 l d T cc& T pS T ] (16) 

C 4 = [ pco T qco T rco T ] (17) 

c 5 = [ 1 e <p s T ] (is) 

where &j in C 5 is an engine throttle parameter. 


These basis functions are designed to model the unknown nonlinearity that exists in the un- 
known plant model. For example, the aerodynamic force in the x- axis for an aircraft can be 
expressed as 


F, - 


firTmax + 2 pV^S (Cl 0 + Cg a 0i + Cg^p + Cg^ui + Ci g S^j Oi 

- 2 * ° v2s ( c Do + C D„« + C-DfiP + + C Ds 5^J (19) 


where the engine thrust is replaced by 8jT max and T max is the maximum engine thrust. 

Thus, Ci, C 2 , and C 3 are designed to model the product terms of a, /3, lo, and 8 in the aero- 
dynamic and propulsive forces. Similarly, C 4 models the cross-coupling terms of the aircraft 
rates in the moment equations, and C 5 models the effects the gravity and propulsive force. 
Alternatively, the basis function Y can also be formed from a subset of C;, i = 1, 2, . . . , 5. 

The update law for the neural net weights W is due to Rysdyk and Calise iRvsd vk & Caliset 
1199 81 and is given by 

W = -T (Ye T PB + ]i ||e T PB|| w) (20) 

where T = r T > 0 _€ R mxm is an adaptive gain matrix, }i > 0 6 R is an e-modification 
parameter ^Narendra & Annaswamvl.ll987 ). . || is a Frobenius norm, and P = P T > 0 € R 6x6 
solves the Lyapunov equation 

PA + A t P = — Q (21) 

for some positive-definite matrix Q = Q T > 0 € R 6x6 . 

The goal is to compute Fi, P 2 , and G by a model inversion indirect adaptive law. The in- 
direct adaptive law updates the estimates of F\, F 2 , and G so that the model inversion con- 
troller 8 can accommodate as much as possible the effects of the unknown plant dynamics. 

Two approaches are considered: 1) an indirect adaptive law based on the Lyapunov's direct 
method, and 2 ) a recursive least-squares indirect adaptive law for parameter estimation of the 
unknown plant model. Both of these approaches are described as follows: 

2.1 Lyapunov-Based Indirect Adaptive Law 

The hybrid adaptive control with model inversion adaptation can be implemented by the 
following indirect adaptive law 

4> = -A(©e T PB + >/||e T PB||o) (22) 

where<J> T = [ WjJ Wj Wj ] £ R 3x k is a weight matrix, © = [ w T Y ( J cr T Yj ] T £ 

RP xl is an input matrix of state and control vectors, A = diag (T w , T a , T^) >06 Rf' x P is an 
adaptive gain matrix, and 1 j = diag >06 Rk x P is an e-modification parame- 

ter matrix. 


Then the estimates of Pi, P 2 , and G can be computed as 


Fi = F* + Wj Y w 

(23) 

F 2 = FZ + WjY a 

(24) 

G = G* + WjY s 

(25) 


The basis functions Y w , Y cr, and Y^ are designed to model the nonlinearity in the plant model 
error. For example, if the plant model error is given by 


e = Aico + A 2 K.L 0 + A^ficv 


( 26 ) 


thenWj = [ A 1 A 2 A 3 ] and Y a , = [ I cd pi ] T . 

The tracking error then becomes 

e = Ae + BW t T + B<1> T © - Be (27) 


The indirect adaptive law {22) can be shown to provide a stable estimation of the unknown 
plant matrices F 1 , P 2/ and G as follows: 

Proof: The matrix A is Hurwitz. Let W = W* + W and <I> = O* + <3> where the asterisk 
symbol denotes the ideal weight matrices that cancel out the unknown plant model error e 
and the tilde symbol denotes the weight deviations. The ideal weight matrices are unknown 
but they may be assumed constant and are bounded to stay within a A-neighborhood of the 
plant model error e, assuming that the input or the command 5 C e Goo is bounded. Then 

A = sup |w* T Y + <l>* T ©-e| (28) 

( o,(t,8 

Choose the following Lyapunov candidate function 

V = e T Pe + tr^W T T“ 1 W + a> T A“ 1 4>) (29) 

where tr (.) denotes the trace operation. 

The time derivative of the Lyapunov candidate function is computed as 

V = e T Pe + e T Pe + 2tr + O t T (30) 

which upon substitution yields 

V = e T (PA + A t p) e + 2 e T PB (w T Y + 0 T © - e) 

+ 2tr [-W T (Ye^PB + fi ||<’ T Pb|| w) - <I> t (®e T PB + i ] ||e T PB|| (31) 

Utilizing the trace operation tr (XY) = YX, where X is a column vector and Y is a row vector, 
then 

2tr (-W T Ye T PB) 

2tr (^-4> T ©e T PB) 

Completing the square yields 

2tr ||e T PB || (W* + W)] = -2 /( ||e T f 

< —}i || e T P B 1 1 (||W|| 2 -||W*|| 2 ) (34) 


= — 2e T PBW T Y (32) 

= — 2e T PB<I> T © (33) 



2 


2tr - 


[-<I> T J7 |e T PB|| (<D* + 4>)] < -2 ||e 

— ^ max (jj') 


pb ; 


c£>* 


^■min (ty) | ^ 2 ^ ^ 

<-||e T PB| [A min (^)||<I>|| 2 -A m „(r ? )||0*|! 2 ] (35) 


where ||.|| is a Frobenius norm, and A m ,„ and A maj - are the maximum and minimum eigenval- 
ues, respectively. 

Then, substituting back into V gives 

V < -e T Qe + 2e T PBA-/Je T PB 11 / ' I|T "H 2 " TA7 *" 2 


) 

- ||e T P B 1 1 [A min (r/) \m 2 -\max(v) ||®*|| 2 ] (36) 


Since ||B|| = 1, it can be established that 

V < -A min ( Q ) if + ||P( fll [2 || A|| - p (||W|| 2 - ||W*|| 2 ) 

-*min(ri)\m 2 + *max(v)\\<S>*\\ 2 ] (37) 

which can also be expressed as 

v < - Hi {a min (Q) |M| - ||P|| [2 || A || + p ||w*|| 2 + A max (r,) ||<l>* || 2 ] 

+ f( ||P||||W|| 2 + A mm (i||P||||<I.|| 2 } (38) 


Let 5 be a compact set defined as 

S = {(e, W, <h) : A mi „ (Q) ||e|| + p ||P|| || W|| 2 + \ min (r,) ||P|| ||<I»|| 2 < r} (39) 

where 

r=||P|| [2 || A || + j.i || TV* || 2 + /\,nax (tj) II c E > * l| 2 ] (40) 

Then V < 0 outside the compact set S. Also there exist functions <pi, y >2 6 K.1Z where 

‘Pi (Hi , II W|| , ||$||) = A min (P) ||c|| 2 + A mi „ (r- 1 ) ||w|| 2 + \ min (a- 1 ) ||<h|| 2 (41) 

P >2 (Hi / ||kV|| / Ill’ll) = A mflx (P) ||e || 2 + A max (r _1 ) || W|| 2 + A max (a -1 ) ||<I>|| 2 (42) 

such that 

<Pi (IklMI^IMI^II) <V<V> 2 (Ili JW|| / ||4>||) (43) 


Then, according to Theorem 3.4.3 of lloannu & Sunl.fl99fih . the solution is uniformly ultimately 
bounded. Therefore, the hybrid adaptive control results in stable and bounded tracking error; 

i.e., e, W,<3? € £co- 


It should be noted that the bounds on ||e||, 1 1 W 1 1 , and ||4>|| depends on ||A||. To improve 
the tracking performance, the magnitudes of A must be kept small. This is predicated upon 
how well the neural network can approximate the nonlinear uncertainty in the plant dynam- 
ics. Increasing the adaptive gains T and A improves the tracking performance but at the same 
time degrades stability robustness. On the other hand, the values of /( and ij must also be kept 
sufficiently large to ensure stability robustness, but large values of ft and t / can degrade the 
tracking performance. Thus, there exists a trade-off between performance and robustness in 
selecting the adaptive gains T and A and the e-modification parameters ft and 7/. 

To ensure that the indirect adaptive law will result in a convergence of the estimates Fj , F 2 , 
and G to their steady state values, the input signals must be sufficiently rich to excite all fre- 
quencies^! i nterest in the plant dynamics. This condition is known as a persistent excitation 
(PE) (ioannu & Sunil 19961) . 

2.2 Recursive Least-Squares Indirect Adaptive Law 

The tracking error equation (TlTb can be expressed as 

e = Ae + + B ^<h T 0 — ej (44) 

Suppose the plant model error can be written as 

e = d> - cv* + Ae = <t> T 0 (45) 

where Ae is the estimation error of A tv. Then, the estimated plant model error is 

e = d>-u>* = a-F^co-F^a-G*S (46) 

where a) is the estimated acceleration. 

The model inversion adaptation using the recursive least-squares indirect adaptive law is 
given by 

4> = -\r& [ e T - 0 T O>) (47) 

m A V / 

R = — \rROQ T R (48) 

m z 

where R = R T > 0 € R px l l is a positive definite covariance matrix and m 2 is a normalization 
factor 

m 2 = 1 + 0 t R© (49) 

The recursive least-squares indirect adaptive law can be derived as follows: 

The estimation error can be minimized by considering the following cost function 


( 50 ) 


To minimize the cost function, the gradient of the cost function with respect to the weight 
matrix is computed and set to zero, thus resulting in 


VjJ 



- 0 T <h) dr = 0 


Equation {51) is then written as 


1 

m 2 


L 


1 0© T dT$ = 
o m z 


f 

Jo 


©e T dr 


(51) 


(52) 


Let 


R- 1 


1 ft T 

/ 00 T dr > 0 
m A Jo 


Differentiating Eq. {53l yields 

dR- 

dt 

It is noted that 

R~ 1 R = I =4 



dR - 1 

dt 


r + r-'r = 0 


(53) 


(54) 


(55) 


Solving for R yields Eq. {48). 

Also, differentiating Eq. {52) yields 

R- 1 ® + -4© T <1> = ^©e T (56) 

m A m A 

Solving for 4> yields the recursive least-squares indirect adaptive law {4 7) . 

The recursive least-squares indirect adaptive law can be shown to provide a stable estima- 
tion of the unknown plant matrices Ti, F 2 , and G as follows: 

Proof: The steady state ideal weight matrix <t>* is assumed to be bounded by a A<j>-neighborhood 
where 

A = sup |<h* T ©-e| (57) 

co,a,8 

The ideal weight matrix W* is assumed to be bounded inside a neighborhood where 

A = sup I W* t Y + 0* T © - i - Ael < sup I W* T Y - Ael + A (58) 

u),(r,8 co,a,8 


Choose the following Lyapunov candidate function 

L = e T Pe + tr ^W t T _1 W + 


( 59 ) 


The only difference between L and V is in the last term. Then, the time rate of change of the 
Lyapunov candidate function is computed as 


L = 


-e T Qe + 2 e T PB (w Tl F + <D T © - e - Ae) 

[w T ( x ¥e T PB + ;; |e T PB|| w)] 

4rd> r © (e T -0 T 0) +0 1 

V / 


— 2tr 


+ tr 


dR- 

~df 


-O 


(60) 


Further simplification yields 


L < -e T Qe + 2e T PBA + 2e T PB$ T @ + p ||e T J PB|| (||W*|| 2 



- ^0 T 4><f> T © + -4 (e T - 0 T 

m z m A \ 



(61) 


L is then bounded by 

L< -Kin ( Q ) lkl! 2 + ll^ll lk|| (2||A||+2 ||<J> t ©|| +/< ||W*|| 2 ) 

- P Ill’ll IMI ||W|| 2 - i ||<I. T ©|| 2 + A ||0 T 0|| ||A|| (62) 

which can also be expressed as 

L<- Ikll [a min (Q) IMI - \\P\l (2 II A || + 2 |0 T ©|| + ft ||W*|| 2 ) 

+P IMII II W|| 2 ] - ^ ||<*> T ©|| (||<*> T 0|| - 2 ||A||) (63) 


L < 0 if 
and 

^ min (q) ikii + nil’ll 


|0 T ©|| >2||A|| (64) 

> ||P|| (2||A||+2||<I> t 0||+ / (||W*|| 2 ) 

> Ill’ll (2 II A || +4 ||A|| + }i ||W*|| 2 ) (65) 


Let C be a compact set defined as 

C = {{e,W,®):A min (Q)\\e\\+ F \\P\\\\W\\ 2 <ror ||o T ©|| < 2 ||A||} 


where 

r=||P|| (2 || A|| +4 ||A|| + }i || W*|| 2 ) 


( 66 ) 

(67) 


Then L < 0 outside the compact set C, and so according to Theorem 3.4.3 of l lloannu & Sunl, 
Il996h . the solution is uniformly ultimately bounded. Therefore, the hybrid adaptive control 
results in stable and bounded tracking error; i.e., e, W, 4> 6 Coo- Thus, the recursive least- 
squares indirect adaptive law is stable. 

The parameter convergence of the recursive least- squares depends on the persistent excita- 
tion condition on the input signals dloannu & Sunl.ll996h . The update law for the covariance 
matrix R has a very similar form to the Kalman filter with Eq. (48} as the differential Riccati 
equation for a zero-order plant dynamics. The recursive least-squares indirect adaptive law 
can also be implemented in a discrete time form with various modifications such as with an 
adaptive directional forgetting factor dBobal et all 1200 5) according to 


%U = <*>; + -4~Ri+ 1©/ W + 1 - ©, T *il (68) 

m f+ 1 L J 

R i+ i = Ri - (w+i + &+ 1) 1 Ri®i®jRi (69) 

where i p and £ are defined as 

fi+l = «?+i - 1 (70) 

to+i = fi + i - 5T 1 (i - ft+ 1) (7i) 


The directional forgetting factor cp is calculated as 


<Pi +- 1 


1 + (1 + p) In (1 + £i+ 1 ) + 


1i+l (1 + fy+l) _ ^ 

. 1 + £i+ 1 + J/i+i 


£i+ 1 
1 + £i+l 


(72) 


where p is a constant, and f/ and are parameters with the following update laws 


1 II T II 2 

Vi+1 = A ;+1 |e;+l - 3>i ©; | 

(73) 

% l = Vi+ 1 (1 + $i) 

(74) 

r il t ii 2 ’ 

^Jt+l = <Pi+ 1 Ajt + (1 + £;+ 1) |£;+1 — < J > ; ©i | 

(75) 


3. Flight Control Simulations 

3.1 Generic Transport Model 

To evaluate the hybrid adaptive flight control method, a simulation was conducted using 
a NASA generic transport mod el (G TM) which represents a notional twin-engine transport 
aircraft as shown in Fig. 2 dlordan et al.l.l200j) . An aerodynamic model of the damaged aircraft 
is created using a vortex lattice method to estimate aerodynamic coefficients and derivatives. 
A damage scenario is modeled corresponding to a 28% loss of the left wing. The damage 
causes an estimated C.G. shift mostly along the pitch axis with Ay = 0.0388c and an estimated 
mass loss of 1.2%. The principal moment of inertia about the roll axis is reduced by 12%, while 
changes in the inertia values in the other two axes are not as significant. Since the damaged 
aircraft is asymmetric, the inertia tensor has all six non-zero elements. This means that all the 
three roll, pitch, and yaw axes are coupled together throughout the flight envelope. 



Fig. 2. Left Wing Damaged Generic Transport Model 


A level flight condition of Mach 0.6 at 4572 m is selected. Upon damage, the aircraft is re- 
trimmed with T = 0.0705W, a = 5.9°, £ = -3.2°, S a = 27.3°, S e = -0.5°, 6 r = -1.3°. The 
remaining right aileron is the only roll control effector available. In practice, some aircraft can 
control a roll motion with spoilers, which are not modeled in this study. The reference model is 
specified by oj v = 2.3 rad/ sec, (v q = 1.7 rad/ sec, cv r = 1.3 rad/ sec, and = £ r = 1/ \/2. 
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The state space model of the damaged aircraft is given by 
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Fig. 3. Pitch Rate 
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Fig. 4. Roll Rate 


The pilot pitch rate command is simulated with a series of ramp input longitudinal stick com- 
mand doublets, corresponding to the reference pitch angle between —3.1° and 3.1°. The track- 
ing performance of the baseline flight control with no adaptation versus the three adaptive 
control methods is compared in Figs. 3 to 6. With no adaptation, there is a significant over- 
shoot in the ability for the baseline flight control system to follow the reference pitch rate as 
shown in Fig. 3. The performance progressively improves first with the direct adaptive con- 
trol alone, then with the hybrid Lyapunov-based indirect adaptive control, and finally with 
the hybrid recursive least-squares (RLS) indirect adaptive control. The Lyapunov-based indi- 
rect adaptive control performs better than the direct adaptive control alone as expected, since 
the presence of the Lyapunov-based indirect adaptive law further enhances the ability for the 
flight control system to adapt to damage. 
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Fig. 5. Yaw Rate 
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Fig. 6. Tracking Error Norm 


The most drastic improvement is provided by the hybrid RLS indirect adaptive control which 
results in a very good tracking performance in all three control axes. In the pitch axis, the 
hybrid RLS indirect adaptive control tracks the reference pitch rate very accurately. In the roll 
and yaw axes, the roll and yaw rate responses are maintained close to zero. In contrast, both 
the direct adaptive control and the hybrid Lyapunov-based indirect adaptive control improve 
the roll and yaw rate responses, but the response amplitudes are still significant and therefore 
can be objectionable particularly in the roll rate. 

Figure 6 is the plot of the tracking error norm for all the three angular rates to demonstrate 
the effectiveness of the hybrid adaptive control method. The hybrid Lyapunov-based indirect 
adaptive control reduces the tracking error by roughly half of that with the direct adaptive 


control alone and by a factor of three when there is no adaptation. Moreover, the hybrid 
RLS indirect adaptive control drastically reduces the tracking error by more than an order of 
magnitude over those with the direct adaptive control and with the baseline flight control. 
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Fig. 7. Bank Angle 




Fig. 8. Angle of Attack 


The attitude responses of the damaged aircraft are shown in Fig. 7 to 9. When there is no 
adaptation, the damaged aircraft exhibits a rather severe roll behavior with the bank angle 
ranging from almost —40° to 20°. The direct adaptive control improves the situation signifi- 
cantly and cuts down the bank angle to a range between about —30° and 10°. With the hybrid 
RLS indirect adaptive control, the bank angle is essentially maintained at its trim value. 


The angle of attack as shown in Fig. 8 is in a reasonable range. The angle of attack when 
there is no adaptation goes through a large swing from 1° to 9°, but the hybrid RLS indirect 
adaptive control reduces the angle of attack to a range between 3° and 8°. 

Figure 9 shows the plot of the sideslip angle. In general, flying with sideslip angle is not a 
recommended practice since a large sideslip angle can cause an increase in drag and more 
importantly a decrease in the yaw damping. With no adaptation, the largest negative sideslip 
angle is about —3°. This is still within a reasonable limit, but the swing from —3° to 1° can 
cause objectionable handling qualities. With the hybrid RLS indirect adaptive control, the 
sideslip angle is retained virtually at zero. 



Fig. 9. Sideslip Angle 


The control surface deflections are plotted in Figs. 10 to 12. Because of the wing damage, 
the damaged aircraft has to be trimmed with a rather large aileron deflection. This causes 
the roll control authority to severely decrease. Any pitch maneuver can potentially run into a 
control saturation in the roll axis due to the pitch-roll coupling that exists in a wing damage 
scenario. With the maximum aileron deflection at 35°, it can be seen clearly that a roll control 
saturation is present in all cases, being the worst when there is no adaptation and the best 
with the hybrid RLS indirect adaptive control. The range of aileron deflection when there 
is no adaptation is quite large. As the aileron deflection hits the maximum position limit, it 
tends to over-compensate in the down swing because of the large pitch rate error produced by 
the control saturation. Both the direct adaptive control alone and the hybrid Lyapunov-based 
indirect adaptive control alleviate the situation somewhat but the control saturation is still 
present. The hybrid RLS indirect adaptive control is apparently very effective in dealing with 
the control saturation problem. As can be seen, it results in only a small amount of control 
saturation, and the aileron deflection does not vary widely. The hybrid RLS indirect adaptive 
control essentially enables the aileron to operate almost at its full authority, whereas with the 
other control methods, only partial control authority is possible. 





Fig. 10. Aileron Deflection 


Figure 11 is the plot of the elevator deflection that shows similar elevator deflections to be 
within a range of few degrees for all the four different controllers. This implies that the roll 
control contributes mostly to the response of the damaged aircraft. 






Fig. 11. Elevator Deflection 


The rudder deflection is shown in Fig. 12. With no adaptation, the rudder deflection is quite 
active, going from —5° to 0°. While this appears small, it should be compared relative to the 
rudder position limit, which is usually reduced as the airspeed and altitude increase. The 
absolute rudder position limit is ±10° but in practice the actual rudder position limit may be 
less. Therefore, it is usually desired to keep the rudder deflection as small as possible. The 
direct adaptive control results in a maximum negative rudder deflection of —4° and the hybrid 
Lyapunov-based indirect adaptive control further reduces it to —2°. The hybrid RLS indirect 


adaptive control produces the smallest rudder deflection and keeps it to less than ±0.5° from 
the trim value. 
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Fig. 12. Rudder Deflection 


3.2 Piloted Flight Simulator 

The Crew- Vehicle System Research Facility (CVSRF) at NASA Ames Research Center houses 
two motion-based flight simulators, the Advanced Concept Flight Simulator (ACFS) and the 
Boeing 747-400 Flight Simulator for use in human factor and flight simulation research. The 
ACFS has a highly customizable flight simulation environment that can be used to simulate a 
wide variety of transport-type aircraft. The ACFS employs advanced fly-by-wire digital flight 
control systems with modern features that can be found in today's modern aircraft. The flight 
deck includes head-up displays, a customizable flight management system, and modern flight 
instruments and electronics. Pilot inputs are provided by a side stick for controlling aircraft 
in pitch and roll axes. 

Recently, a piloted st udy has be en conducted in the ACFS to evaluate a number of adap- 
tive control methods dCampbell et alJ , l2010l) . A high-fidelity flight dynamic model was devel- 
oped to simulate a medium-range generic transport aircraft. The model includes aerodynamic 
models of various aerodynamic surfaces including flaps, slats, and other control surfaces. The 
aerodynamic database is based on Reynolds number corrected wind tunnel data obtained 
from wind tunnel testing of a sub-scale generic transport model. The ground model with 
landing gears as well as ground effect aerodynamic model are also included. 

A number of failure and damage emulations were implemented including asymmetric dam- 
age to the left horizontal tail and elevator, flight control faults emulated by scaling the control 
sensitivity matrix (B-matrix failures), and combined failures. Eight different NASA test pilots 
were requested to participate in the study. For each failure emulation, each pilot was asked 
to provide Cooper-Harper Ratings (CHR) for a series of flight tasks, which included large 
amplitude attitude capture tasks and cross-wind approach and landing tasks. 



Fig. 13. Advanced Concept Flight Simulator at NASA Ames 



Fig. 14. Pilot Evaluation of Adaptive Flight Control 




Seven adaptive control methods were selected for the pilote d study that include e-modification 
fNarendra & Ann aswamy|. 1987) , hybr id adaptive control iNguven et alJ.l2006h . optimal con- 


trol modification (‘Ngugeneta, 


stability method (NgugenetaL 


2008|), metric-driven adap tive control using bounded linear 


j 20071), C\ adaptive control ICao & H ovakimvanl.l200Sl). adap- 
tive loop recovery iCalise et alJ.l2009t) . and composite adaptive control iLavretsky[ 20091) . This 
is by no means an exhaustive list of new advanced adaptive control methods that have been 
developed in the past few years, but this list provides an initial set of adaptive control meth- 
ods that could be implemented under an existing NASA partnership with the industry and 
academia sponsored by the NASA Integrated Resilient Aircraft Control (IRAC) project. 


The study generally confirms that adaptive control can clearly provide significant benefits 
to improve aircraft flight control performance in adverse flight conditions. The study also 
provides an insight of the role of pilot interactions with adaptive flight control systems. It 
was observed that many favorable pilot ratings were associated with those adaptive control 
methods that provide a measure of predictability, which is an important attribute of a flight 
control system design. Predictability can be viewed as a measure of how linear the aircraft 
response is to a pilot input. Being a nonlinear control method, some adaptive control methods 
can adversely affect linear behaviors of a flight control system more than others. Thus, while 
these adaptive control methods may appear to work well in a non-piloted simulation, they 
may present potential issues with pilot interactions in a realistic piloted flight environment. 
Thus, understanding pilot interaction issues is an important consideration in future research 
of adaptive flight control. 


Cooper-Harper Rating Differentials 



Fig. 15. Cooper-Harper Rating Improvement of Various Adaptive Control Methods 


With respect to pilot handling qualities, among the seven adaptive flight controllers evaluated 
in the study, the optimal control modification, the adaptive loop reco very, and th e composit e 
adaptive control appeared to perform well over all flight conditions dCampbell et al.i.i2010l) . 
The hybrid adaptive control also performs reasonably well in most cases. For example, with 
the B-matrix failure emulation, the average CHR was 5 for 8 capture tasks with the baseline 
dynamic inversion flight controller. The average CHR number was improved to 3 with the 




hybrid adaptive control. In only one type of failure emulations that involved cross-coupling 
effects in aircraft dynamics, the performance of the hybrid adaptive flight controller fell below 
that for the e-modification which is used as the benchmark for comparison. 

Future NASA research in advancing adaptive flight control will include flight testing of some 
of the new promising adaptive control methods. Previously, NASA conducted flight test- 
ing of the Intelhgent^hght i Contro^ < (IFQ ii o^ < a i NASA < FA5 i aircraft < u£^ntil i ^008j[Bosworth 
& Williams-Hayes, 20071) . In' January of 2011, NASA has successfully completed a flight test 
program on a NASA F-18 aircraft to eva lua te a new adaptive flight controller based on the 
optimal control modification ^Nguyen et alj.|200 8!). Initial flight test results indicated that the 
adaptive controller was effective in improving aircraftSs performance in simulated in-flight 
failures. Flight testing can reveal new observations and potential issues with adaptive control 
in various stages of the design implementation that could not be observed in flight simulation 
environments. Flight testing therefore is a critical part of validating any new technology such 
as adaptive control that will allow such a technology to transition into production systems in 
the future. 

4. Conclusions 

This study presents a hybrid adaptive flight control method that blends both direct and in- 
direct adaptive control within a model inversion flight control architecture. Two indirect 
adaptive laws are presented: 1) a Lyapunov-based indirect adaptive law, and 2) a recursive 
least-squares indirect adaptive law. The indirect adaptive laws perform on-line parameter es- 
timation and update the model inversion flight controller to reduce the tracking error. A direct 
adaptive control is incorporated within the feedback loop to correct for any residual tracking 
error. 

A simulation study is conducted with a NASA wing-damaged transport aircraft model. The 
results of the simulation demonstrate that in general the hybrid adaptive control offers a po- 
tentially promising technique for flight control by allowing both direct and indirect adaptive 
control to operate cooperatively to enhance the performance of a flight control system. In 
particular, the hybrid adaptive control with the recursive least-squares indirect adaptive law 
is shown to be highly effective in controlling a damaged aircraft. Simulation results show 
that the hybrid adaptive control with the recursive least-squares indirect adaptive law is able 
to regulate the roll motion due to a pitch-roll coupling to maintain a nearly wing-level flight 
during a pitch maneuver. 

The issue of roll control saturation is encountered due to a significant reduction in the roll 
control authority as a result of the wing damage. The direct adaptive control and the hybrid 
adaptive control with the Lyapunov-based indirect adaptive law restore a partial roll control 
authority from the control saturation. On the other hand, the hybrid adaptive control with the 
recursive least-squares indirect adaptive law restores the roll control authority almost fully. 
Thus, the hybrid adaptive control with the recursive least-squares indirect adaptive law can 
demonstrate its effectiveness in dealing with a control saturation. 

A recent piloted study of various adaptive control methods in the Advanced Concept Flight 
Simulator at NASA Ames Research Center confirmed the effectiveness of adaptive control in 
improving flight safety. The hybrid adaptive control was among the methods evaluated in the 


study. In general, it has been shown to provide an improved flight control performance under 
various types of failure emulations conducted in the piloted study. 

In summary, the hybrid adaptive flight control is a potentially effective adaptive control strat- 
egy that could improve the performance of a flight control system when an aircraft operating 
in adverse events such as with damage and or failures. 
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